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Abstract 

Dyson-Schwinger equations provide a prominent approach to physics of 
strong interactions. To reproduce the hadronic phenomenology well, 
the Dyson-Schwinger approach in the rainbow-ladder approximation 
must employ an effective interaction between quarks which is fairly 
strong at intermediate ( Q 2 ~ 0.5 GeV 2 ) spacelikc transferred mo¬ 
menta. We have recently proposed that such an interaction may orig¬ 
inate from the dimension 2 gluon condensate (A 2 ) which has recently 
attracted much attention, and showed that the resulting effective run¬ 
ning coupling leads to the sufficiently strong dynamical chiral symme¬ 
try breaking and successful phenomenology at least in the light sector 
of pseudoscalar mesons. In the present paper, we give a more detailed 
investigation of the parameter dependence of these results. 

Pacs: 11.10.St; 11.30.Qc; 12.38.Lg; 14.40.Aq 

1 Introduction 

In recent years, the dimension 2 gluon condensate (A“ A a ^) = (A 2 ) attracted 
a lot of theoretical attention [1-7], to quote just several of many papers offer¬ 
ing evidence that this condensate may be important for the nonperturbative 
regime of Yang-Mills theories, particularly QCD. Although (A 2 ) is not gauge 
invariant, it was even argued that its value in the Landau gauge may have 
a physical meaning [2,3,7]. In our recent paper [8] we argued that (A 2 ) 
may be relevant for the Dyson-Schwinger (DS) approach to QCD. Namely, 
in order that this approach leads to a successful hadronic phenomenology, 
an enhancement of the effective quark-gluon interaction seems to be needed 
at intermediate ( Q 2 ~ 0.5 GeV 2 ) momenta 1 , and Ref. [8] showed that the 

x We adopt the convention k 2 = —Q 2 < 0 for spacelike momenta k. 
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gluon condensate (A 2 ) provides such an enhancement. It also showed that 
the resulting effective strong running coupling leads to the sufficiently strong 
dynamical chiral symmetry breaking and successful phenomenology in the 
light sector of pseudoscalar mesons. However, the issue of the parameter 
dependence of the results was just commented on very briefly. Thus, in the 
present paper, in Sec. 03 we give a more detailed investigation and presen¬ 
tation of the parameter dependence of these results. A brief recapitulation 
of the DS approach and the effective interaction it needs is given in the next 
section. 

2 DS approach and its effective interaction 

DS approach to hadrons and their quark-gluon substructure [9-11] has 
strong and clear connections with QCD. Besides being covariant, this ap¬ 
proach is chirally well-behaved and nonperturbative. This has been crucial, 
especially in the light-quark sector of QCD, for successful descriptions of 
bound states achieved by phenomenological DS studies (e.g., see recent re¬ 
views [10,11] and references therein), where one can treat soundly even the 
processes influenced by axial anomaly 2 , which is really remarkable for a 
bound-state approach. What happens is that in the process of solving DS 
equations, one in essence derives a constituent quark model which turns out 
to be successful over a very wide range of masses. Its chief virtue is that it 
incorporates the correct chiral symmetry behavior through the gap equation 
for the full, dynamically dressed quark propagator S q and the Bethe-Salpeter 
(BS) equation for the bound states of the dynamically dressed quarks (and 
antiquarks). That is, the constituent quarks arise through dressing result¬ 
ing from dynamical chiral symmetry breaking (DySB) in the (“gap”) DS 
equation for the full quark propagators, while the light qq pseudoscalar 
solutions of the BS equation (in a consistent approximation) are (almost 
massless) guasi-Goldstone bosons of DySB. Generation of DySB is well- 
understood [9,10,18-22] in the rainbow-ladder approximation (RLA). Thus, 
phenomenological DS studies have mostly been relying on RLA and using 
Ansatze of the form 

[K(k)} h e ‘ = iiimM-k 2 ) D‘ b „(k) 0 [L YU [y VW (1) 

for interactions between quarks. In this equation, e, /, g, h schematically rep¬ 
resent spinor, color and flavor indices and D^ u {k )o is the free gluon propaga- 

2 See, e.g., Refs. [12,13] for the 7r° —> 77 transition amplitude Tff, and Refs. [14-17] 
for the related transition 7 —> 7r + 7r°7r~. 
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Figure 1: The effective non-strange (q = u ) quark mass function M u (—Q 2 ) 
calculated using the effective coupling ® proposed in Ref. [8] and the input 
parameters given by Eqs. m and PI). 

tor in the gauge in which the aforementioned DS studies have been carried 
out almost exclusively, namely the Landau gauge: 

= ^(-!V + . (2) 

while a e s(Q 2 ) is an effective running coupling on which we will comment 
below at length. 

The BS equation for the bound-state vertex T^/ of the meson composed 
of the quark of the flavor q and antiquark of the flavor q', is then 

/ d^P P P 

The consistent RLA requires that the same interaction kernel m be previ¬ 
ously used in the DS equation for the full quark propagator S q . That is, 
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( 4 ) 


dressed quark propagators S q (k) for various flavors q, 

Sq 1 ^) = A q{p 2 )i> - B q(p 2 ) , (<? = U , d , S, ...) , 

are obtained by solving the gap DS equation 

/ d^f \ a \ b 

aeff(-(p -t) 2 )D“ b u (p -£)0 Y YS q {£)— Y . 

( 5 ) 

Following the approach of Munczek and Jain [19,20], the gap equation © 
is unrenormalized, but regularized by an ultra-violet cutoff L. This cutoff 
is however huge compared to the QCD scale A qcd- (In the present paper, 
L = 134 GeV as in Ref. [20].) In Eq. ©, rn q is the cutoff-dependent bare 
mass of the quark flavor q breaking the chiral symmetry explicitly. The case 
fh q = 0 corresponds to the chiral limit where the current quark mass m q = 0, 
and where the constituent quark mass M q ( 0) = B q (0)/A q (0) stems exclu¬ 
sively from the nonperturbative phenomenon of D\SB. Of course, calling 
the “constituent mass” the value of the “momentum-dependent constituent 
mass function” M q (p 2 ) = B q (p 2 )/A q (p 2 ) at exactly p 2 = 0 and not at some 
other low — p 2 , is a matter of a somewhat arbitrary choice. However, it is 
just a matter of terminology and nothing essential. What is important to 
get a successful hadronic phenomenology, especially in the light-quark sector 
(q = u,d,s), is that DySB is sufficiently strong. This means that the gap 
equation © should yield quark propagator solutions A q (p 2 ) and B q {p 2 ) giv¬ 
ing the dressed-quark mass function M q (p 2 ) whose values at low —p 2 are of 
the order of typical constituent mass values, namely several hundred MeV, 
even in the chiral limit. A typical example of such M q (p 2 ) is given in Fig. 
12 obtained with a e s(Q 2 ) © proposed originally in our Ref. [8] and further 
advocated in the present paper. 

Indeed, the issue of the origin of the interaction 0. or, equivalently, 
a e s(Q 2 ) which would enable successful phenomenology is crucial for the DS 
studies. The form of a e ff is only partially known from the fact that at large 
spacelike momenta it must reduce to a pe rt (Q 2 ), the well-known running 
coupling of perturbative QCD. However, for momenta Q 2 < 1 GeV 2 , where 
non-perturbative QCD applies, the interactions are still not known; there¬ 
fore, in phenomenological DS studies, a e s(Q 2 ) must be modeled for Q 2 < 1 
GeV 2 - e.g., see Refs. [9-11,20-23]. There, one can see that phenomeno¬ 
logically most successful of those modeled interactions have a rather large 
bump at the intermediate momenta, around Q 2 ~ 0.5 GeV 2 . For example, 
in Fig. [2 compare a e s(Q 2 ) used by Jain and Munczek (JM) [20] and by 
Maris, Roberts and Tandy (MRT) [10,11,21,22]. In any case, successful 
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DS phenomenology requires that this modeled part of the interaction 0 be 
fairly strong. That is, regardless of details of the interaction, its integrated 
strength in the infrared must be fairly high to achieve acceptable description 
of hadrons, notably mass spectra and DySB [10,11]. 

Theoretical explanations on what could be the origin of so strong nonper- 
turbative part of the phenomenologically required interaction are obviously 
very much needed, either from the ah initio studies of sets of DS equations 
for Green’s functions of QCD (see, e.g., the recent review [9]) or from some¬ 
where outside DS approach. The particularly important result of the ab 
initio DS studies is that, in the Landau gauge, the effects of ghosts are ab¬ 
solutely crucial for the intermediate-momenta enhancement of the effective 
quark-gluon interaction [9,24-27]. This is obvious in the expression for the 
strong running coupling a s (Q 2 ) in these Landau-gauge studies [9,24-27], 

a B (Q 2 ) = a s (n 2 )Z(Q 2 )G(Q 2 ) 2 , (6) 

where a s (fa 2 ) = g 2 /4 tt and Z(fa 2 )G(fi 2 ) 2 = 1 at the renormalization point 
Q 2 = /a 2 . The gluon renormalization function Z(—k 2 ) defines the full gluon 
propagator D^ b u (k) in the Landau gauge: 

D;i(k) = Z(-k 2 )D;i(k )o = 6 ab [~ 9lxv + . (7) 

Similarly, G(—k 2 ) is the ghost renormalization function which defines the 
full ghost propagator (/c) = S ab G(—k 2 )/k 2 . 

While the ab initio DS studies [9, 24-27] do find significant enhance¬ 
ment of a s (Q 2 ), Eq. 0, until recently this seemed still not enough to yield 
a sufficiently strong DySB (e.g., see Sec. 5.3 in Ref. [9]) and a success¬ 
ful phenomenology. However, for carefully constructed dressed quark-gluon 
vertex Ansatze, Fischer and Alkofer [24] have recently managed to obtain 
good results for dynamically generated constituent quark masses and pion 
decay constant f n , although not simultaneously also for the chiral quark- 
antiquark ( qq } condensate, which then came out somewhat larger than the 
phenomenological value. Thus, the overall situation is that there is progress 
in this direction [24-28], but that further investigation and elucidation of 
the origin of phenomenologically successful effective interaction kernels re¬ 
mains one of primary challenges in contemporary DS studies [10,11]. This 
provided the motivation for our paper [8], where we pointed out that such an 
interaction kernel for DS studies in the Landau gauge resulted from cross¬ 
fertilization of the DS ideas on the running coupling of the form 0 [9,24-27] 
and the ideas on the possible relevance of the dimension 2 gluon condensate 
(A^AW) = (A 2 ) [1-7,29-32], 


5 





Figure 2: The momentum dependence of various strong running couplings 
mentioned in the text. JM [20] and MRT [10, 22] a e s(Q 2 ) are depicted 
by, respectively, dashed and dash-dotted curves. The effective coupling © 
proposed and analyzed in the present paper is depicted by the solid curve, 
and a s (Q 2 ) © of Fischer and Alkofer [24] (their fit A) by the dotted curve. 
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In Ref. [8], we gave arguments that the (H 2 ^contributions to the OPE- 
improved gluon (H) and ghost (G) polarization functions (found a long time 
ago by Refs. [29-32] and more recently confirmed by Kondo [4]) lead to an 
effective coupling a e s(Q 2 ) given by 

«eff(Q 2 ) = « pert (Q 2 ) Z NpeTt (Q 2 ) G Npert (Q 2 ) 2 , (8) 


where a pert (Q 2 ) is the running coupling of perturbative QCD, and 


Z Npert (Q 2 ) 


G Npert (Q 2 ) 


1 

i i i c A 
1 ^ W ^ W 

1 


1 - 


m A , Co 

w w 


( 9 ) 

( 10 ) 


The functions Z Npert (Q 2 ) and G Npert (Q 2 ) are the nonperturbative (Npert) 
parts of the, respectively, gluon and ghost renormalization functions Z(Q 2 ) 
and G(Q 2 ). They crucially depend on the quantity iua which can be in¬ 
terpreted as a dynamically generated effective gluon mass, and which is 
proportional to the dimension 2 gluon condensate (A 2 ). Concretely, for the 
Landau gauge (to which we stick throughout this paper), the number of 
QCD colors N c = 3 and the number of space-time dimensions D = 4, 

m 2 A = g 2 {A 2 ) = -m 2 c , ( 11 ) 

where me is a dynamically generated effective ghost mass. (In a subsequent 
work, Kondo et al. [5] also worked out logarithmic corrections to Eq. (1111) 
thanks to which the dynamical gluon mass (and ghost mass) vanishes as 
Q 2 —» oo, as it must according to, e.g., Cornwall [33,34], However, taking 
this into account is not necessary at the degree of refinement and precision 
at which we work in this paper.) 

For g 2 (A 2 ), the Landau-gauge lattice studies of Boucaud et al. [1] yield 
the value 2.76 GeV 2 . This is compatible with the bound resulting from the 
discussions of Gubarev et al. [2,3] on the physical meaning of (A 2 ) (although 
it is gauge-variant) and its possible importance for confinement. We thus 
use this value in Eq. m and obtain 


mA = 0.845 GeV . 


( 12 ) 


In our considerations below, this value will turn out to be a remarkably good 
initial estimate for the dynamical masses rriA and mg. 
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The coefficients Ca and Cg appearing in Z Npert (Q 2 ) Q and G Npert ((5 2 ) 
m, can, in principle, be related to various other condensates [30-32], 
but some of them are completely unknown at present. Therefore, both 
Ca and Cg should at this point be treated as free parameters to be fixed 
by phenomenology. Fortunately, Ref. [8] managed to make the estimate 
Ca = (0.640 GeV) 4 . This estimate [8] is based on the role of only one 
condensate [35], the well-known gauge-invariant dimension 4 condensate 
(F 2 ) [36], and thus misses some (unknown) three- and four-gluon contri¬ 
butions [31,32], Therefore, and since the true value of (F 2 ) is still rather 
uncertain [37], we do not attach too much importance to the above precise 
value of Ca but just use it as an inspired initial estimate. 

There is no similar estimate for Cg, but one may suppose that it would 
not differ from Ca by orders of magnitude. We thus try 

C G = C A = (0.640 GeV) 4 (13) 

as an initial guess. It turns out a posteriori that this value of Cg leads to a 
very good fit to phenomenology. 

As we discussed in Ref. [8], Eq. ® can be justified for relatively high 
Q 2 , but not for low Q 2 . For example, a pe rt (Q 2 ) must ultimately hit the 
Landau pole as Q 2 gets lowered. However, this can be handled as in 
other phenomenological DS studies. Their various choices of a e s(Q 2 ) usu¬ 
ally also contain a pei -t(Q 2 ), but since handling the Landau pole problem at 
the fundamental level is out of their scope, they [20-22, 38-40] just shift 3 
the Landau pole to the timelike momenta in all logarithms appearing here: 
hi(Q 2 /Aq CD ) —> ln(xo + Q 2 /A.q CD ). Presently, we adopt this latter proce¬ 
dure. Concretely, for a pert (Q 2 ) we use throughout the MS-scheme two-loop 
expression used before by JM [20] and our earlier phenomenological DS 
studies [38-42], This means we use throughout the infrared (IR) regulator 
xq = 10 (to which all results are almost totally insensitive), the number 
of quark flavors Nf = 5, and A qcd = 0.228 GeV. These parameters of 
expert {Q 2 ) are thereby fixed and do not belong among variable parameters 
such as Ca, Cg, the variation of which is discussed below. 

In the present context, the more serious objection to our a e g (JHJ) is that 
we cannot in advance give an argument that the factor Z Npert (Q 2 ) G Npert (Q 2 ) 2 
in the proposed a e s(Q 2 ) (JHJ) indeed approximates well nonperturbative con¬ 
tributions at low Q 2 (say, Q 2 < 1 GeV 2 ), but can only hope that our results 

3 As pointed out already by, e.g., Cornwall [33], dynamically generated gluon mass can 

provide the physical reason for such a change in the arguments of logarithms. That is, 
x 0 oc m 2 A /A 2 QCD ~ 10. 
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to be calculated will provide an a posteriori justification for using it as low 
as Q 2 ~ 0.3 GeV 2 [since Eq. ® takes appreciable values down to about 
Q 2 ~ 0.3 GeV 2 ]. Of course, Z Npert (Q 2 ) and G Npert (<5 2 ) must be wrong in 
the limit Q 2 —> 0, as they are based on the results derived by OPE [4,29-32], 
which certainly fail in that limit. For example, detailed investigations of the 
Q 2 — * 0 asymptotic behavior in ab initio DS studies [9,24-27], settled down 
to the conclusion that a s (Q 2 ) remains finite as Q 2 —> 0, which is also sup¬ 
ported by several lattice calculations [43,44]. On the other hand, if the 
presently interesting (A 2 ) condensate is explained by an instanton liquid, the 
coupling vanishes as a s (Q 2 ) oc Q 4 [6], which is closer to the behavior of our 
a e s(Q 2 ) ©• Still, Eqs. HU enforce, for small Q 2 , even much more dramatic 
suppression of our a e g(Q 2 ) (JHJ), which vanishes as Q 12 . This is an unrealistic 
artefact of the proposed form (JHJ) when applied down to the Q 2 —> 0 limit. 
Nevertheless, because of the integration measure in the integral equations 
in DS calculations, integrands at these small Q 2 [where our a e g(Q 2 ) ® 
is doubtlessly too suppressed] do not contribute much, at least not to the 
quantities (such as (qq) condensate, meson masses, decay constants and am¬ 
plitudes) calculated in phenomenological DS analyses. Hence, the form of 
a es{Q 2 ) Q 2 close to zero is not very important 4 for the outcome of these 
phenomenological DS calculations. This is because the most important for 
the success of phenomenological DS calculations seems the enhancement at 
somewhat higher values of Q 2 - e.g., see the humps at Q 2 ~ 0.4 to 0.6 GeV 2 
in the JM [20] or MRT [21, 22] a e g(Q 2 ), dashed curves and dash-dotted 
curves in Fig. [2j Our a e g(Q 2 ) ® exhibits such an enhancement centered 
around Q 2 ~ m\/ 2, as shown by the solid curve representing it in Fig. |21 
This enhancement is readily understood when one notices that Eq. ® has 
four poles in the complex Q 2 plane, given by 

(Q 2 ) 1,2 = ^ m 2 A =F i\J 4 C G - m\ ^ [poles of G Npert (<5 2 )] , (14) 

(Q 2 ) 3,4 = ^ ^ ~rni\ T i\J 4 Ca -m^ [poles of Z Npert (Q 2 )] . (15) 

For min-fOG, Ca} > rn\ /4 there is no pole on the real axis, but a saddle 
point in the middle of two complex conjugated poles. For the DS studies, 
which are almost exclusively carried out in Euclidean space, spacelike k 2 (i.e., 
Q 2 > 0 in our convention) is the relevant domain and is thus pictured in Fig. 
121 There, the maximum of a e g{Q 2 ) (JHJ) at the real axis is at Q 2 ~ m 2 A / 2, 

4 Of course, the Q 2 — > 0 domain would give an important contribution in a case with a 

sufficiently strong (but still integrable) divergence in a e ff(Q 2 ), such as the delta function 

in Ref. [23]. 
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i.e., the real part of its double poles (<3 2 )i, 2- The height and the width of 
the peak is influenced by both Cq and m A - The enhancement of a e s(Q 2 ) 
© is thus crucially determined by the (A 2 ) condensate through Eq. (1111) . 
and by the manner this condensate contributes to the ghost renormalization 
function, which enters squared into the effective coupling ©• 

3 Phenomenology with the condensate-enhanced 
coupling 

We solved the DS equations for quark propagators and BS equations for 
pseudoscalar qq (q = u, d , s) bound states in the same way as in our previous 
phenomenological DS studies [38-41]. This essentially means as in the JM 
approach [20], except that instead of JM’s a e s(Q 2 ), Eq. © is employed 
in the RLA interaction ©. We can thus immediately present the results 
because we can refer to Refs. [38-41] for all calculational details, such as 
procedures for solving DS and BS equations, all model details, as well as 
expressions for inputs such as the aforementioned IR-regularized a per t(Q 2 ) 
and explicit expressions for calculated quantities, e.g., for f n . 

3.1 In the chiral limit 

In the chiral limit, where the bare (and current) quark masses vanish, the 
only parameters are those defining our a e s(Q 2 ) ©, namely tua, Ca and Ce¬ 
lt turns out that the initial estimates m and m, motivated above, need 
only a slight modification to provide a very good description of the light 
pseudoscalar sector: it is enough to increase the estimate tua = 0.845 GeV 
by just 5%. That is, the parameter set 

C A = (0.640 GeV) 4 = C G , rn A = 0.884 GeV (16) 

leads to (to begin with) an excellent description of DySB, which gives rise 
to Goldstone bosons which are also massless pseudoscalar qq bound states. 
This is seen in the first line of Table© our good chiral limit values of the pion 
decay constant ( f n & 88 MeV) and the qq condensate [(qq) ~ (—214 GeV) 3 ] 
satisfy the Gell-Mann-Oakes-Renner (GMOR) relation (two last columns 
in Table© very well, at the level of a couple of percent. These chiral-limit 
results are similar to, e.g., the corresponding results with JM a e s(Q 2 ), which 
are also given in Table © (in the last line) for comparison. 

The behavior of the momentum-dependent constituent mass function 
M q (p 2 ) = B q (p 2 ) / A q (p 2 ) is also qualitatively similar both to M q (p 2 ) found 
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«eff, Cq, 

C a , rn A 

(qq) [GeV 3 ] 

U [GeV] 

_<ff> [GeV] 

J 7r 

lim ™^o S t GeV l 

Eqs. ©,(©2) 

(-0.214) 3 

0.0882 

1.261 

1.293 

Eqs. ©,(HZD 

(-0.217) 3 

0.0905 

1.241 

1.289 

JM CKeff [20] 

(-0.227) 3 

0.0898 

1.368 

1.401 


Table 1: The chiral-limit results for f n and ( qq) and the test of the GMOR 
relation for our a e ff © and the JM one [20]. The quark condensate and the 
current quark mass m are calculated at the renormalization scale /r = 1 GeV. 
In DS approach, good values of f n automatically lead to good description 
of 7T° —> 77, since the empirically successful amplitude = l/ATT 2 f n is 
always obtained analytically in this approach in the chiral limit [12,13]. 


earlier by JM [20] and ourselves [38-42] with JM a e fj and to M q (p 2 ) obtained 
now with our a e g © but with different parameters (this is exemplified by 
M u (-Q 2 ) in Fig. CJ. Quantitatively, for the parameters m and the chiral 
limit (fhq = 0), the constituent quark mass M q ( 0) = 0.306 GeV. This is 
almost 25 % below both our old results for M u ( 0) [39] obtained with JM 
a e g and our present M u ( 0) in Fig. IB pertaining to the refitted parameters 
(11711 and (much less importantly) to m u / 0 (THU) . However, the quantitative 
differences of such a size are not a problem, since calculations in practice 
show that a successful reproduction of the hadronic phenomenology require 
just that values of this (anyway unobservable) quantity at low Q 2 are of the 
order of several hundred MeV, i.e., of the order of typical constituent quark 
mass, M q ( 0) ~ M nucleon /3 ~ M p / 2. 

The constituent quark mass in the chiral limit, directly related to the (qq) 
condensate, is also very convenient for illustrating the dependence of the key 
D\SB phenomenon on the model parameters. If we vary Cq (for fixed values 
of niA and Ca) away from its phenomenologically favorable value in Eq. 
m, which gives sufficient enhancement of a e ff, the dynamically generated 
constituent quark mass M q ( 0) quickly falls. Beyond some critical value of 
Cq, it is always exactly zero, meaning that the D\SB is then completely 
absent. The sensitivity of our results to Cq is understandable, since from 
Eqs. m it is clear that Co, in combination with tua, influences the height 
and width of the peak of a e s(Q 2 ) © for spacelike momenta. In spite of 
this sensitivity, we were able to find other combinations of parameter values 
which lead to good results. For example, the values 

C A = (0.6060 GeV) 4 = C G , m A = 0.8402 GeV (17) 
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Figure 3: The dependence of the dynamically generated constituent quark 
mass M q (0 ) on the parameter Cq illustrates the disappearance of DySB for 
unfavorable values of Cq- when for given values of rriA and fh q (here rriA = 
0.884 GeV and fh q = 0) Cq deviates from the value that gives sufficient 
enhancement of a e ff, the dynamically generated mass M q ( 0) quickly falls. 
Moreover, beyond some critical value of Cq, it is always exactly zero since 
the DySB phenomenon then completely disappears. 

yield the second line of Table Cl This indicates that there may be an in¬ 
teresting interplay between rriA and Cq and motivates us to find how the 
phenomenologically favorable values of rriA and Cq are related. However, 
we will do it below in the more realistic, massive case, away from the chiral 
limit. There, the quark bare masses (and the related current masses) deviate 
from zero so that empirical masses of pseudoscalar mesons can be obtained. 

3.2 Away from the chiral limit 

We start by noting that both of the two sets of (mi, Ca, Cq) values quoted 
above as successful in the chiral limit, Eqs. m and m, gives a good fit 
also away from the chiral limit. As the first shot, we adopt without any 
change the explicit breaking of chiral symmetry from JM, that is, the bare 


12 



fleffi Cg, Ca, 
tua, fh u , m s 

H 

M h [MeV] 

f H [MeV] 

Ty [MeV” 1 ] 

Eqs. ©,<nnj> 

7r 

136.70 

91.2 

0.272 x 10- a 

and (11811 

K+ 

520.72 

112.1 


Eqs. ©.(HI 

7r 

136.17 

93.0 

0.256 x 10 " 3 

and (11811 

K+ 

516.28 

112.5 


Eqs. ©,(T7} 

7r 

134.96 

92.9 

0.256 x 10 " 3 

and mi 

K+ 

494.92 

111.5 


experimental 

7 T U 

134.9766 ± 0.0006 

91.9 ±3.5 

(0.274 ± 0.010) x 10 " 3 

values 

K+ 

493.677 ± 0.016 

112.8 ± 1.0 



Table 2: The masses and decay constants of pions and kaons, and the 7 r° —> 
77 decay amplitude TJ7, obtained in DS approach with our a e fr(Q 2 ) (| 8 |). The 
hrst two lines result from the initial parameters rn a , Ca,g © and the quark 
bare mass parameters m hxed already by the broad JM phenomenological 
fit [ 20 ]. These masses (1TK1) with another (niA, Ca,g) parameter set m give 
the third and the fourth line. Similarly, the fifth and the sixth line result 
from a e s(Q 2 ) with Ca,g given by Eq. inn, and the slightly altered bare 
masses ©■ The last two lines are the corresponding experimental values. 

mass parameters (fh q ) of light quarks (q = u, d, s ) leading to the broad 
phenomenological fit with their a e g- [ 20 ], namely 

fh u = rhd = 3.1 ■ 10 ~ 3 GeV , fh s = 73 • 10 ~ 3 GeV . (18) 

These values of fh U) d lead to an excellent description of the pion as a quasi- 
Goldstone boson of Dy;SB also in conjunction with our a e fj © and Eq. m , 
as witnessed by the hrst line in Table [21 where we predict pion mass, weak 
decay constant, and 7 r° —> 77 amplitude very close to their empirical values 
(in the seventh line of Table ©. For the same reason as in the chiral limit, 
the results are again quite sensitive to changes of Cq but not to Ca- Table 01 
illustrates this relatively weak sensitivity to the changes of Ca for the case of 
the parameter set m&m- Table 01 shows one can increase (or decrease) 
C A by a factor of 2 , and the results change little. 

The second line of Table 01 reveals that the parameter set 
works somewhat less well in the strange sector, as the kaon mass is 5% too 
high. However, a deviation of this size is not worrisome in the present cir¬ 
cumstances where we know that the model interaction anyway misses some 
aspects (such as the Q 2 —> 0 behavior and non-ladder contributions), and 
where we just want to point out that the (A 2 ) condensate is a possible source 
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set A 

set B 

H 

M h 

Ih 

M h 

Ih 

7 r 

K+ 

ss 

0.132609 

0.490013 

0.716287 

0.0877711 

0.107034 

0.131859 

0.136671 

0.499146 

0.727431 

0.0955129 

0.113735 

0.133728 


Table 3: This table illustrates rather weak sensitivity to changes of Ca- The 
“set A” of input parameters is given by Eqs. d and HZD with the change 
Ca ‘ZCa- The “set B” of input parameters is given by Eqs. with the 
change Ca Ca/2- Meson masses and meson decay constants are in units 
of GeV, while ss stands for the non-physical pseudoscalar ss bound state. 

of the needed enhancement of a e g(Q 2 ). In fact, the empirical success in the 
strange sector is quite reasonable considering that we used the standard 
JM mass parameters [20], (as we did also in [38-41]) and no refitting was 
performed there (although a e s(Q 2 ) was different). 

Nevertheless, it is interesting to see what changes are brought by refit¬ 
ting. If one for example tries the values of tua^Cg and Ca given by Eq. 
C2> instead of Eq. (fTH , one gets the third and the fourth line in Table El 
instead of, respectively, the first and second line. Thus, the improvement 
achieved thereby is not significant, indicating that we should try changes of 
the bare quark masses fh q . It turns out that slight changes of the values ca 
are sufficient to achieve agreement with experiment in the both non-strange 
and strange sectors. For example, the parameter set which gives the fifth 
and sixth lines of Table El thus reproducing the empirical mass of both 7 r° 
and K + together with good results for their decay constants and 7 r° —> 77 
amplitude Tj,7, is given by mi, Cq and Ca from Eq. ca and by the bare 
quark masses 

fh u = m d = 3.046 ■ 10 " 3 GeV , fh s = 67.70 • 10 " 3 GeV . (19) 

This parameter set, Eqs. ca and ca, is also the one giving the gap 
equation solutions resulting in the momentum-dependent constituent mass 
function M q (—Q 2 ) displayed in Fig. Cl 

The parameter set ca&ca also gives us a good description of the 77 - 77 ' 
complex, along the lines of our Refs. [39,42], Although it means employ¬ 
ing just a minimal extension of the DS approach, we must relegate this to 
another paper [45]. 

The preferred parameter set ca&ca is a result of a systematic exam¬ 
ination of refitting possibilities performed by studying the dependence on 
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2.8 


2 . 9 


3.1 3.2 3.3 


Figure 4: The curves are the solutions of the equations F = 2.5% and 
F = 5.0% in ( rh u ,rh s ) plane. The point is the position of the simple, non¬ 
degenerate minimum at the bare quark mass values Pl). 


the input parameters x = ( m. u ,m s ,mA,CG,CA ) of the function 



2 


x 100% 


( 20 ) 


namely the sum of squared differences of the four experimentally measured 
(Vexp) and presently theoretically calculated (yth) quantities y G {M n o, f n ±, M K o, _/%±}. 
We kept choosing Ca = Cq for simplicity, since we find that moderate vari¬ 
ations of Ca do not affect our results much anyway, as already illustrated 
by Table [3 

Minimization of Eq. ca shows different respective characters of the 
a e ff parameters (mu, Cq, Ca) and the mass parameters ( fh u , fh s ). The point 
m in the parameter subspace (m u , m s ) is the location of a non-degenerate 
minimum of the function P|) . Thus, the possible values of the bare quark 
masses ( fh u , m s ) can be precisely restricted by demanding that the function 
(f20|) be below certain value. Figure 01 shows F = 5.0% and F = 2.5% curves 
in the (fh u ,m s ) plane, with mu and Cq{= Ca) fixed at Eq. (1T71) . At the 
minimum, for (fh u ,fh s ) values m, we obtain F ss 1.5%. 

In contrast to the bare quark masses (fh u ,fh s ), the parameters defining 
a e fi cannot be determined so unambiguously. By this we do not mean just 
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0.836 0.838 0.84 0.842 0.844 


Figure 5: F vs. (m^C^/ 4 ) contour plot. The darkest color corresponds to 
F ~ 1.5%, defining the valley of the minimal F. Conversely, the lighter the 
shade of gray, the larger the value of F, i.e., the overall difference between 
the calculated and experimental quantities. 
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the aforementioned weak sensitivity to Ca- They also cannot be fixed by 
minimization of F m in the same sense as the bare quark masses even 
though the results are very sensitive to rriA and Cq- The point is that F has 
no simple minimum in the (rriA, Cq^)— plane as it has in ( fh u ,fh s ) plane: Fig. 
0 reveals a minimum in the form of a narrow, straight “valley” described 
very well by a linear relation between m.A and Cq A . Thus, in spite of high 
sensitivity to rriA and C G , there are many pairs of these quantities which 
give a fit comparable (within few percent) to that resulting from the values 
m, as long as they approximately satisfy the linear relation 

(C G ) 1/4 = 0.7742 m A - 0.0444 GeV . (21) 

That is, the function (PH) measuring the difference between the calculated 
and experimental values of M n o ,f n ±, M K o , f K ± has a degenerate minimum 
in the shape of a narrow valley. It is bounded by the values (Cq) m i n ~ (0.6 
GeV) 4 and (Ca) m a Y ~ (0.9 GeV) 4 in the sense that between these values we 
managed to find solutions providing excellent fits (F of the order 1.5%) to 
the empirical values. 

4 Conclusion 

The dimension 2 gluon condensate (A 2 ) enabled the derivation [8] of a suit¬ 
ably enhanced a e g (Q 2 )- This effective interaction leads to the sufficiently 
strong DySB and successful phenomenology at least in the light sector of 
pseudoscalar mesons. This opens the possibility that instead of modeling 
a eS (Q 2 ), its enhancement at intermediate Q 2 may be understood in terms 
of gluon condensates, which seem to provide an important mechanism pro¬ 
posed and studied for the first time in our recent Ref. [8]. The systematic 
examination of the parameter space, i.e., various fitting possibilities set forth 
in the present paper, allows us to conclude that this scenario is compatible 
with reasonable values of both (A 2 ) -condensate and the gauge-invariant di¬ 
mension 4 gluon condensate (F 2 ) [36]. In the relevant momentum region, 
ctes(Q 2 ) (and thus also the solutions of DS and BS equations and results for 
calculated measurable quantities) depend only very weakly on Ca, which 
parametrizes contributions of dimension 4 condensates to the gluon propa¬ 
gator. The essential parameters Cq and m-A, on which the dependence is 
very strong, are not independent. Thus, due to the relation (121(1 . Eq. © 
is an essentially one-parameter model for a e ff, albeit on a relatively small 
interval of Cq- This can be interpreted as another instance that what counts 
is the integrated strength of the interaction. Over the possible range, we 
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have a continuous set of parameter pairs (via-, Cg); their values are such 
that they give higher peaks at smaller squared momenta, resulting in simi¬ 
lar integrated strengths. We find that the phenomenologically allowed range 
of values of the dynamically generated gluon mass via is in agreement with 
the lattice results [1] on ( A 2 ) in the Landau gauge. Also, phenomenologi¬ 
cally allowed values of Cq, which parametrizes contributions of dimension 
4 condensates to the ghost propagator, are such that they might be a sign 
that Cq is indeed mostly determined by the dimension 4 gluon condensate 
(F 2 ) [36], 
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